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1 Introduction 



The standard model (SM) of elementary particles, despite it's spectacular success in 
explaining Physics upto scale close to TeV, suffers from the well known gauge hierarchy 
problem, where the Higgs boson mass receives quantum loop corrections of the order of 
Planck scale. To restrict the mass of the Higgs boson within electroweak scale we have to 
fine tune these quantum corrections at each loop by adjusting the counter terms in the 
SM. Resolution to this problem has been explored in the context of supersymmetric and 
extra dimensional models. 

Another problem encountered within the framework of SM is related to the masses 
of the fermions. The known fermions in the SM have masses ranging from 1 eV to 172 
GeV. The reason why the fermions have different and hierarchical masses is one of the 
challenging problems in particle physics. Various approaches to solve the fermion mass 
hierarchy problem have been discussed over a wide range of works [1]. In this work we 
propose yet another possible solution to this problem in a multiply warped geometry 
model which has already been shown to offer a resolution of the gauge hierarchy problem 
in the context of Higgs mass. 

Higher dimensional models offer a new explanation for stabilizing the Higgs boson mass 
around the electroweak scale. Among these, the model proposed by Randall and Sundrum 
(RS) has a warped geometry with one extra spatial dimension |2]. The extra spatial 
dimension in this model supports two 3-branes and the warping along this dimension can 
generate exponential mass scale suppression ~ 10^^ between these two 3-branes. Hence, 
one of the 3-branes can be interpreted as Planck-brane and the other one to be TeV- 
brane. By confining the Higgs boson on the TeV-brane, we generate TeV scale mass for 
the Higgs field. Subsequently this model was generalized to accommodate more than one 
extra warped dimension and it was qualitatively indicated that the requirements of the 
resolution of gauge hierarchy problem as well as little hierarchical moduli give rise to 
the possibility of resolving the fermion mass hierarchy problem [3]. Here, we demonstrate 
explicitly that in a multiply warped model, the fermion mass hierarchy problem can indeed 
be addressed. 

In the multiply warped braneworld models [3], the number of 3-branes increases with 
the number of extra dimensions. Half of these 3-branes have energy scale clustered around 
TeV scale while the scale of the remaining 3-branes would be around Planck scale. For 
example, in a 6-dimensional doubly warped model, there will we four 3-branes, out of 
which two will have scale around TeV and the remaining two would be around the Planck 
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scale. By computing the fermion Yukawa couplings on different 3-branes clustered around 
TeV scale, we show that the respective masses of fermions will split. Based on this, we 
argue that in a multiply warped braneworld model, the fermion mass splitting between 
these 3-branes can offer an explanation of the observed fermion mass hierarchy of the SM. 
More specifically, we take the 6-dimensional warped model as our case study and we show 
that by adjusting the parameters of the model to appropriate values, the fermion mass 
splitting on the two TeV-branes can be generated. By extending this analysis to more 
than six dimensions, we show the possibility of explaining all fermion masses of the SM. 

The paper is organized as follows. In the next section we give a brief overview of 
multiply warped braneworld models. In Sec. 3 we analyze the fermion mass splitting on 
two different 3-branes in the 6-dimensional warped model. We generalize this to higher 
dimensional multiply warped braneworld models, and argue that the different fermion 
masses in the SM can be explained from the mass splittings that are generated in these 
models. We conclude in Sec. 4. 

2 Multiply warped braneworld model 

In this section we give brief description of six and seven dimensional warped spacetime 
models [3], which are extended from the 5-dimensional RS model. Further generalization 
of these models beyond seven dimensions can be worked out in a similar fashion [3]. 

The six dimensional warped model has the usual four spacetime and two extra spatial 
dimensions. The two extra spatial dimensions are successively compactified on S1/Z2 
orbifold. The manifold of the six dimensional spacetime is [M^'^ x /Z2\ x /Z2- The 
metric ansatz consistent with this manifold is 

ds^ = b\z)[a\y)7]^^dx^'dx'' + Rldy^] + r^dz^, (1) 

where x^, /i = 0, ■ ■ ■ , 3, are the four non-compact spacetime coordinates, y and z are the 
angular coordinates of the two extra dimensions with the moduli Ry and r^, respectively. 
a{y) and b{z) are some functions of the extra coordinates. Here, rj^i, = diag(— 1, 1, 1, 1). 
Since the y and z coordinates are orbifolded on S^/Z2, with the orbifold fixed points 
y = 0,TT and z = 0,7t, there are four 4-branes at these fixed points. Intersection of any 
of the two 4-branes give one 3-brane. Hence, in this model there would be four different 
3-branes located at {y, z) = (tt, it), (tt, 0), (0, vr) and (0, 0). We take one of these 3-branes 
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as our Tev scale universe. The total bulk-brane action of this model has a form [3] 

S = 5*6 + 5*5, 5*6 = y" d'^xdydzy/-g6{Re - A), 

S5 = [ d^xdydz[Vi5{y) + V25{y - tt)] + ! d^xdydz[V^5{z) + Vi5{z - tt)] (2) 



Here, Vi,2 and V3,4 are brane tensions of the branes located at ?/ = 0, tt and 2; = 0, tt, 
respectively. A is the cosmological constant in 6-dimensions. 

After solving Einstein's equations with orbifolded boundary conditions the solution to the 
warp functions of the metric of Eq. ([1]) has a form [3] 

cosh(A;z) 



a{y) = exp{-c\y\), b{z) 



cosh{kTT) 



Ryk , / -A 



^ " r,cosh(fc7r)' ^"""^VlOMl,- 

Here, Mp is the Planck scale. 

The functional form of a{y) in Eq. ([3]) is similar to the exponential warping in the 
5-dimensional RS model |2]. In addition the function b{z) gives warping along z direction. 
It can be noticed that for the 3-brane located at y = 0, z = tt we have no warping and 
hence it is called Planck brane while for the 3-brane located at y = vr, 2; = the warping 
is maximum and is called standard model or TeV brane. However, from the relation 
for c in Eq. (|3]), in order not to have large hierarchy between Ry and r^, one of c or 
k must be large and other is small. This indicates large warping in one direction and 
small along the other. It is now easy to show that when Eq. ([3]) is considered together 
with the maximum warping condition at the visible brane i.e. a(7r)6(0) = 10^^® then the 
minimum hierarchy between the two moduli can be achieved for c ~ 11 and k ~ 0.4. Now, 
with this choice of values for c and k, the 3-branes located at {y, z) = (vr, 0) and (tt, tt) get 
exponential suppression ~ 10^^^ with respect to the 3-branes located at {y, z) = (0, 0) and 
(0,7r). However due to small warping along z-direction, the mass scales on the 3-branes 
at {y, z) = (0, 0) and (0, vr) would of the order of Planck scale while the mass scales on the 
3-branes at {y, z) = (vr, 0) and (vr, vr) would be suppressed close to TeV scale. So in the 
six dimensional warped model, the requirement of absence of any large hierarchy between 
the two moduli along with both the moduli implies that the scales of two 3-branes are 
clustered around TeV scale and that of the other two are at around Planck scale. The 
moduli of such six dimensional models can be stabilized to their desired values following 
the mechanism suggested by Goldberger and Wise where a bulk scalar field is introduced 
as a stabilizing field with appropriate vacuum expectation values at the boundaries [1] 
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We now describe the seven dimensional model which has three extra spatial dimen- 
sions. Each of the extra dimension is orbifolded by symmetry. The manifold under 
consideration of this model is [{M^'^ x IZ2\ x 5*^/^2] x /Z2. Denoting the three extra 
dimensions by angular coordinates z and w, the bulk-brane action in this model can 
be written as 

S = Sf + Sq, Si = J d'^xdydzdw^J—g-j{R-j — A7), 

Sq = J d'^xdydzdw[Vi6{w) + ¥26(10 — vr)] + J d'^xdydzdw[V36{z) + Vi6{z — vr)] 

+ j d^xdydzdw[V55{y) + Ve5{y-7r)]. (4) 

Here, A7 is cosmological constant in 7-dimensions. Vi^...^ are brane tensions of 5-branes 
located at if = 0, vr, z = 0, vr and y = 0,tt, respectively. The metric solution of the above 
action which satisfies Einstein's equations and satisfies the orbifolded boundary condition 
is [3J 

= ( '''''I'alf''! [exp{-2cy)r]^,dx>'dx'' + Rldy^] + r^dz^] + Rldw^ 

cosh [In) { cosh (ktt) J 



15 ' Rw cosh(Z7r) 



(5) 



R^ cosh(A;7r) cosh(/7r) cosh(A;7r) 

Here, R^ is the moduli corresponding to the w direction. Looking at the above solution of 
metric, we have exponential suppression along the y direction, whereas, the suppression 
along z and w are given by cosine-hyperbolic functions. Once again a consistent set of 
values for the dimensionless quantities would be: c ~ 10, A; ~ 0.1 and / ~ 0.1, so that 
there will not be large hierarchies among the moduli Ry, and R^,. 

In the seven dimensional model, the 5-branes located at w = 0, vr, 2; = 0, vr and y = 0,tt 
intersect and give twelve different 4-branes which lie at {z, w) = (0, 0), (0, vr), (vr, 0), (vr, vr), 
{z, y) = (0, 0), (0, vr), (vr, 0), (vr, vr) and {y, w) = (0, 0), (0, vr), (vr, 0), (vr, vr). Intersection of 4- 
branes results in 3-branes, and in this model there are eight different possibilities, which lie 
at {y,z,w) = (0,0,0), (0,0, vr), (0,vr,0), (vr, 0,0), (vr,vr,0), (vr,0,vr), (0,vr,vr) and (vr,vr,vr). 
For the consistent set of values of c, k and I, which are described in the previous paragraph, 
the 3-branes located at (0, 0, 0), (0, 0, vr), (0, vr, 0) and (0, vr, vr) have no suppression, hence, 
these 3-branes can be identified as Planck scale branes. Whereas, the 3-branes located at 
(vr, 0, 0), (vr, 0, vr), (vr, vr, 0) and (vr, vr, vr) have exponential suppression ~ 10^^ to TeV scale 
masses, and hence they may be identified as TeV scale branes. So in the seven dimensional 
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model, four 3-branes are clustered around TeV scale and the remaining four are clustered 
around Planck scale. 

The six and seven dimensional braneworld models, which are described here, can easily 
be generalized to higher dimensions. A notable feature in all these models is that half 
of the total number of 3-branes would be clustered around TeV scale and the remaining 
3-branes would be clustered around Planck scale. As explained in Sec. 1, we exploit 
this feature and try to see how much splitting between fermion masses can be obtained 
geometrically. In the next section we describe the calculation of fermion mass splitting 
in the six dimensional model. We can generalize this calculation to higher dimensions, 
and we argue that in a multiply braneworld model the hierarchical mass pattern of SM 
fermions can be explained. 



3 Fermion mass splitting in a multiply warped braneworld 
model 

As explained previously, we take 6-dimensional warped model as our case study and cal- 
culate fermion mass splitting between two different 3-branes. A viable way of generating 
fermion masses is through Yukawa couplings. We generalize the approach taken in [S] , for 
generating fermion masses in multiply braneworld models. We assume that both the Higgs 
boson and left-handed doublet quark fields of SM are confined to the usual 4-dimensional 
world. We allow the right-handed fields of the SM to be in the bulk of 6-dimensions. Con- 
sider the Yukawa couplings for left-handed doublets with its right-handed counter-parts. 
Since right-handed fermions are allowed in the bulk, we can estimate the gauge invariant 
Yukawa couplings on both the 3-branes at (y, z) = (vr, 0), (vr, tt), which are given below. 

5'vis = j d^Xy/-g^is%LoLix)Ho{x)ipR{x, vr, 0) + h.c, 

'S'vis' = J d'^Xy/-g^is'YQLoLix)Ho{x)tpR{x, vr, vr) + h.c. (6) 

Here, Yq, Yg' are Yukawa couplings in 6-dimensions in the two branes. Note that their 
dimensions are M~^, where M ~ Mp. Lql indicates left-handed SU(2) doublet of either 
quark or lepton field of the SM, and ipR is its right-handed counter-part. Hq is the Higgs 
doublet. (?vis = det(5'|^j,(7r, 0)) and g^is' = det(5'j^^(7r, vr)), where gl,j j^ is the metric function 
defined in Eq. ([T]). Since the fermions are 8-component spinors in 6-dimensions, we have 
to choose the representation for Higgs field Hq accordingly in the above equation. Also 



6 



since ipn is a bulk field, it decomposes into Kaluza-Klein (KK) fields in the 4-dimensional 
world. After the Higgs acquires vacuum expectation value, the Yukawa couplings in above 
equation give mixing masses between the left-handed fermion and the KK tower of right- 
handed fields. The zero mode of the KK tower is identified as the SM field. Nevertheless, 
to be consistent we have to include all possible mixings with the KK fermions, and after 
diagonalizing it we can identify the lowest eigenvalue to be the SM fermion mass. In fact, 
this is the procedure employed in [5j, and here we generalize this to 6-dimensions. For 
this we have to study the KK excitations of a bulk fermion field, which is done in for a 
6-dimensional model. However, the nonzero KK masses and their wave functions are not 
computed in [B] . In the next subsection we calculate them which will be necessary for the 
computation of SM fermion masses. 



3.1 KK fermions 

The kinetic part of the action for a bulk fermion field \I' in 6-dimensions is O [71 18] 



Sf = j d^xdydzV^E^ 
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(7) 



Here, the capital letter A denotes index in the curved space and the lower case letters 
a, b, c denote indices in the tangent space. UbcA is spin connection and E^ is inverse 
vielbein. Here, we do not write the bare mass of bulk fermion since the SM fermions 
acquire mass through the Higgs mechanism via Eq. (E]). The spin connection term in the 
above equation give zero contribution due to the diagonal form of the metric in Eq. ([T]). 
The order of the Dirac matrices in 6-dimensions would be 8x8, and they can be taken 
as [9]: 

Yi^ = Y^a'^, r^ = n5^a\ r5 = z75®a' (8) 

Here, are the Dirac matrices in 4-dimensions and 75 = 27°7^7^7^. a*, i = 1,2, 3, 
are the Pauli matrices and cr° is the 2x2 unit matrix. The chirality in 6-dimensions is 
defined by the matrix f = r0rir2p3p4p5 g^Q]^ tlia.t r\l/± = ±\l/±. The chiral fermions 
in 6-dimensions have both left- and right-handed chirality of 4-dimensions, which can be 
projected by the operators Pl,r = (1 T ^r°r^r^r'^)/2. Since, fermions have chirality in 
six dimensions, we can write ip = ip+ + Each field carrying either plus or minus 
chirahty can further be decomposed into left- {4'±l) and right-chiral {i'±R) fields in four 
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dimensions. These fields can be expressed in KK expansion as [6] 

Here, fields of the form ^^^'''\x^) are the KK fields in 4-dimensions and the KK wave 
functions should satisfy the following normalization condition and eigenvalue equations 

dydzh\z)a\y) {ftiU-R,-Liy^ f^/£L,-R,-Liy, ^) = ^^'^'S''"', (10) 

itVy + V,)f^^£\y,z) = -M,,,f^^i\y,z), 

[-rVy + V^)f^lt\y,z) = -M,,,f^^'^\y,z), 

{iVy + V^)f^lt\y,z) = M,,,/^]?)(y,z), 

(-tVy + V,)/_^i\y,z) = M,,k&iy,z), (11) 

where the differential operators: Vy = ^^{4:dya + 2ady), = i^a{bd-J) + 2hdz) and M,- /,. 
is the mass of the KK fermions ^^^•'^K In [B] the KK wave function for zero mode has been 
solved, which is identified with Mo,o = 0. By writing f^j^\y, z) = fy{y)fz{z), we get the 
solution for fy and fz as 

/I \ exp ( ^^£ptan"^(tanh(A;z/2)) cosh(A;7r) ) 

Here, ci is a separation constant which is defined in [6], and can be worked out from the 
normalization of wave function f^^\y, z), using Eq. flTI]]) . The nonzero KK masses and 
their wave functions are necessary to calculate the mixing masses of Eq. (JH]), which will 
be worked out below. 

The eigenvalue equations in Eq. ( ITTj) can be transformed into second order differential 
equations as 



Vy + V,){-^Vy + V,)f<^ff^, = Uy^^il,. (13) 

Here we will show how the nonzero eigenvalues ^ can be found for the wave function 
f+R ■ The wave functions of other chiral modes can be analogously worked out. By writing 
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= x^''''^Ky)^^'^'^K^)^ the above second order differential equation can be recast into 
1 [24c2 - 20cdy + 4^2] x^^-'^) ^ Mfj^ 



where the first and second hnes of the above equation depends on functions of y and 2, 
respectively. Hence, we can choose 

1 [24c2 - 20c9^ + 492]x(^''=) 



1 (59^6 + 269^) ^^^^''^^ 2ic (5(9^6 + 269^)^ 



(15) 



where is a separation constant of mass-square dimensions. By making the transfor- 
mations 

_ g|cy = ^ cosh(A;7r)e"^ k' = k/r„ (16) 

the first of Eq. (ITSi) can be put into a Bessel's equation as 

Zf-^ + + {Zl - .^)x^-'^ = 0, (17) 



where 



2 25 / 2 6c \ cosh (/ctt) 



Now, the solution to the wave function x'^^''^\y) can be written as 

x^'^'\y) = ^[Uz,) + hY^{z,)i (19) 

where and h are some constants. can be determined from the normalization condition 
of Eq. ffTOj) . h and the KK masses Mj^k can be determined by applying the following 
boundary conditions dyX^-'''^\y)\y=o,-K = 0. The boundary condition at ?/ = gives us that 
the constant h is negligibly small and the condition at y = tt gives the following relation 

hUy) + y[J,-i{y) - J.+i(y)] = 0, (20) 

where y = cosh(fc7r)e^'^. The above equation determines the nonzero mass eigenvalues 
Mj of the KK modes for a specific value of u. The index u depends on the separation 
constant m^. This value can be determined by solving the equation along the z, which is 
the second of Eq. (|T5i) . 
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The second of Eq. f|T5l) can be expressed as 

(59,6 + 2hd,)[{bd,h + 269,)e(^'') + = ^2^^^''^ (21) 

where ki = 2icrz/Ry and A;2 = —'irlrn?^. The solution to the above equation can be found 
from the hnear equation {hdzh + 2bdz)^^-''^^ = k^i^^'^\ where fcs satisfies the quadratic 
equation fc3(A;3 + ki) = k2. After doing this, the general solution to the second of Eq. f|T5|) 
can be written as 



_,_ exp(fc^ tan~-'^[tanx.^ 2 



exp(fc^tan ^[tanh( y)] cosh(A;7r)/A;) 

'/^{kz) 
icr. - / ARlm^, . 



cos 



where ai,2 are constants, and one of them can be determined by the boundary condition 
dz$,^^'''\z)\z=o,TT = 0. The condition at z = gives us Oi = —a2k^/k^. And the condition 
at 2; = TT gives the following relation 

KiH - 5A;tanh(A;7r))^+(7r) - kj{k^ - 5A;tanh(A;7r))^"(7r) = 0. (23) 

The undetermined constant 02 can be found from the normalization condition of Eq. 
(fTOj) . The left-hand side of Eq. fl23l) has a sinusoidal behavior, and zeros of this function 
determine the nonzero values of m^, which should be put in Eq. (1201) to find the nonzero 
KK masses Mj^k of fermions. 

3.2 Fermion mass splitting 

After calculating the KK masses and their wave functions we now return to Eq. ([6j), 
and calculate the SM fermion masses on both the 3-branes at {y,z) = (7r,0), (7r,7r). To 
make the kinetic term canonical we need to do rescaling of the SM doublet Lql and the 
Higgs field Hq, since we have confined them in the 4-dimensional world. Imagine that our 
universe is identified as {y, z) = {it, 0), then the kinetic energy terms of a fermion ipo and 
scalar 0o on this brane are 

S = J d^x./^, [gZd4>*od,(^o + E^Hn, 0)V'or'^a^^o + ■ ■ ■ ] (24) 

The inverse vielbein is given by Ej^^rr, 0) = 5^e'^'^ cosh(fc7r) and we have g'^^^ = e'^^'^ cosh^(fc7r)?7'' 
Substituting them in the above equation, we get 



S = d^x 



p — 2c7r p— 3c7r 



cosh {kn) cosh (kTc) 

10 



(25) 



Hence, to get the canonical kinetic terms for scalar and fermion fields in 4-dimension, we 
have to do the following rescalings 



(po = e^'^ cosh(A;7r)</), ipQ 



e 2 cos 



h3/2(^7r)V', 



(26) 



where and if) are renormalized scalar and fermion fields. Recall Sec. 2 where we pointed 
out that c ~ 10 and ~ 0.1 give a consistent set of values. For this set of values 
e'^'^ cosh(A;7r) ~ 10^^, which naturally gives (0) ~ 100 GeV. This is in similar lines of what 
it is shown in the 5-dimensional RS model [5]. 

The SU(2) doublet fermion of the SM can be put in either of the + or the — chirality 
of the 6-dimensional world. Without any loss of generality, we choose it to be in the — 
chirality of 6-dimensional spinor. Then, after substituting the KK wave function expan- 
sion for \E'/j, Eq. and also doing the rescaling of the 4-dimensional fields Lql and Hq, 
the first of Eq. ^ turns out to be 



e 2-- 



X- 



COS 



h=^/2(^^) 



(27) 



Here, = Yg/ yRy^z- Yg is the effective 4-dimensional Yukawa coupling which is an 0{1) 
quantity. After the Higgs field acquires vacuum expectation value (vev), the left-handed 
SM fermion gets mixing masses with all the right-handed KK fields, which can be written 
as 

— — en 

[moh = eosh^/2(fe^/ +^'^(^' 0)^6^;, (28) 

where v is the vev of Higgs field. In the previous subsection, we have seen that the fields 
i/j^'^^ have Dirac masses with their left-handed counter parts. Hence, in the basis "^l = 



of the SM doublet and ipR 



(^i?, '^+r\ '^+r\^+r\ ■■■), where ipi is a component 
- "^^'ri we get the following form for mixing matrix 





^ ["^d]o,o 




["^d]o,i 











Mi,o 










M = 








Mo,i 


















Mi,i ■ 





v 



(29) 



We have numerically computed the nonzero KK masses Mij, whose values are coming 
in the multi-TeV scale. From the relations given in previous subsection, the KK masses 
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depend on apart from k and c. For r^^ ~ 10^^ GeV, the lowest KK mass has been found 
to be around 6 TeV. Here we have fixed k = 0.7 and c = 11.24 where the justification of 
these values will be given later (see Table [1]). Some of the modes like Mo,i, Mo,2, etc. may 
decouple away since their masses and wave functions are found out to be close to zero. 
We have also found that the off-diagonal values of [^dIjj' are coming in the electroweak 
scale. As a result due to the high values of the KK masses, after diagonalizing the above 
matrix, we can see that the zeroth mode mass [m£)]o,o is negligibly affected. Hence, the SM 
fermion mass is approximately the same as the zeroth mode mass [m£)]o.o- The same kind 
of feature can be seen if we repeat this calculation on the other 3-brane at {y, z) = (vr, vr), 
i.e. on the second term of Eq. ([H]). Because of this, to calculate the fermion mass splitting 
between the two 3-branes, {y, z) = (vr, 0), (vr, vr), we work out the corresponding values of 
the zeroth mode. 

From Eq. (!28l) we can notice that the following factor 

" :/IT(-.o) = ^W- (30) 



COS h=^/2(A;7r)'+" ' ' ' cosh3/2(A;7r) 

determines the zeroth mode mass on the 3-brane at [y, z) = {it, 0) where Yqv ~ 100 GeV. 
In the above equation we have used ff^\y,z) = fy{y)fz{z) and also Eq. f[T^ . The 
corresponding factor on the 3-brane at {y, z) = (tt, tt) would be 

F' = e-t'="cosh^/'(fc7r)|/|°^°)(7r,7r)| = e^(^i+^)\ (31) 

The ratio of masses confined on 3-branes at {y, z) = (tt, 0), (tt, tt), is 

^ ^ ^ cosh'/\kn) ^^^^ 

This ratio R gives us the splitting between the fermions on these two 3-branes. The 
mass term of a fermion confined on {y, z) = (vr, 0) brane and the mass term of a fermion 
confined on the other TeV 3-brane at {y, z) = (vr, vr) will have a relative suppression by a 
factor given by R. In Table [H we give possible values of R in the 6-dimensional model. 
Here we have chosen c so that e^'^ cosh(/c7r) = 10^^, which yields the desired suppression 
from Planck to TeV scale in the 6-dimensional model [6]. Table 1 has been constructed 
for different choices for the two moduli in a six dimensional model where the Planck to 
TeV scale warping is achieved for every choice. Such a hierarchical pattern of masses can 
simultaneously be obtained by including larger number of extra dimensions and adjusting 
the different moduli according to Table 1. It has been explained in section 2 that for seven 
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k 


c 


Ry_ 


R 


0.7 
1.2 
1.7 
2.2 
2.7 


11.24 

10.75 

10.24 

9.75 

9.25 


73.3 
194.34 
628.9 
2223.7 
8268.9 


1 

9.75 
1 

101.1 
1 

1065.7 
1 

11243.4 
1 

r\j =- 

10^ 


4.2 


7.75 


495757.0 




6.1 


5.85 


~ 10^ 





Table 1: Various mass splitting ratios R for different values of k and c. ^ gives the 
hierarchy between the two moduli. 

or higher dimensional models the stack of 3-branes have energy scales clustered around 
TeV scale and Planck scale. 

Here, using Table 1, we give few examples of how the fermion mass hierarchy can 
be addressed in multiply braneworld models. Without loss of generality we can localize 
top quark to the TeV-brane at {y, z) = (tt, 0), which has a mass of ^172 GeV. Then the 
fermion mass at the other TeV-brane, i.e. at {y,z) = [it, it), depends on the values of k 
and c. For k = 0.7 and c = 11.24, the fermion field localized at {y, z) = {tt, tt) would have 
a mass of ~10 GeV, which is a reasonable value to explain the bottom quark mass of ;^4.2 
GeV. On the other hand, if we set k = 1.2 and c = 10.75, the fermion mass on the 3-brane 
at {y,z) = (tt, tt) would be suppressed to ~1 GeV, which is in the right ball park region 
to explain either the charm quark mass (~1.3 GeV) or the tau lepton mass (~1.77 GeV). 
Likewise, if we set k = Q.l and c = 5.85, the suppression on the 3-brane at {y, z) = {tt, tt) 
is such that we can even explain the neutrino mass scale of 1 eV. However, the values in 
Table 1 are generated for the case of six dimensional warped model, where we could only 
explain the masses of two fermion fields. By generalizing this to seven dimensions, we will 
have four 3-branes at around TeV scale and we also get additional freedom in the choice 
of the moduli Rw Hence, in the seven dimensional warped model, adjusting the value 
of the other modulus according to Table 1 we can explain four fermion masses. These 
arguments can now easily be generalized in a braneworld model having more than seven 
dimensions where SM fermions of different masses are localized on different 3-branes of 
the stack of clustered 3-branes around Tev scale through the usual localization mechanism 
[10] . The hierarchy in standard model fermion masses thus can be explained. 
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4 Conclusions 



In this work we have proposed a possible origin of the fermion mass hierarchy in the 
standard model through multiple warped geometry. The requirement of the resolution 
of gauge hierarchy problem restricts the energy scale of half of the 3-branes in such 
model around TeV while the remaining half are close to Planck scale. Using the mass 
generation mechanism of the standard model fermions proposed in where the right 
handed fermions are allowed to propagate in the bulk and both the Higgs boson and 
left-handed fermion fields of SM are confined to the 3-branes, we have shown that we can 
adjust and tune various moduli to have different mass splittings once one mass parameter 
is determined by fixing the value of the only unknown parameter Yqd- Assuming that the 
standard model fermions are localized at different 3-branes closely clustered around TeV 
the slight difference in warping results into different contributions from the right handed 
fermion wave function in the Yukawa coupling which in turn leads to different masses of 
the SM fermion. The hierarchical mass pattern of the standard model fermions thus can 
be generated by this mechanism. It is interesting to note from Table 1 that larger is the 
ratio between two moduli more is the splitting of fermion masses which are localized in 
those two branes. The multiple warped model thus proposes a geometric understanding 
of the splitting of the masses of standard model fermions. 
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